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I. INTRODUCTION 

Since its discovery in 2003 the X(3872) charmo- 
nium is subject of many experimental and theoretical ef- 
forts aimed at disclosing its nature — for a recent review 
see [2| . The most recent data on the mass of the X is [H 

Mx = (3871.85 ± 0.27(stat) ± 0.19(syst)) MeV, (1) 

with a width of Fx < 1.2 MeV. However, the problem 
of the quantum numbers for the X is not fully resolved 
yet: while the analysis of the tt+tt" J/?/; decay mode of 
the X(3872) yields either 1++ or 2 ^ quantum numbers 
0,13: the recent analysis of the 7r+7r~7r° J/-0 mode seems 
to favour the 2 ^ assignment Q, though the 1++ option 
is not excluded. This question is clearly very central, 
for the most promising explanations for the X in the S- 
wave DD* molecule model Q as well as in the coupled- 
channel model require the quantum numbers l"*"^. In 
addition, the X cannot be a naive cc 2 ^ state, for its 
large branching fraction for the D^D^tt^ mode [Ml is not 
compatible with the quark- model estimates for the 2 ^ 
charmonium So, for the 2 quantum numbers, very 
exotic explanations for the X would have to be invoked. 

The aim of the present paper is to perform a combined 
analysis of the data on the tt+tt^ and vr+vr^Tr*' mass dis- 
tribution in the 7r+7r~J/'0 and ir^ J / ip mode, re- 
spectively. We find that the S-wave amplitudes from the 
decay of a 1++ state provide a better overall description 
of the data than the P-wave ones from the 2 '", espe- 
cially when the parameter range is restricted to realistic 
values. We conclude then that the existing data favour 
1^"*" quantum numbers of the X, however, improved data 
in the tt+tt^vt'^ J/?/; mode are necessary to allow for defi- 
nite conclusions regarding the X quantum numbers. 



II. EXPERIMENTAL SITUATION 

Recently Belle announced Q the updated results of the 
measurements for the reaction X(3872) tt+tt^ J/V': 

B}^ = [8.61 ± 0.82(stat) ± 0.52(syst)] x 10-^ 



"27r 



[4.3 ± 1.2(stat) ± 0.4(syst)] x 10" 



for the charged J/ippK^ and neutral B'^ — > 

J/ippK'^ mode, respectively, with B2Tr being the product 
branching fraction Br(B -> i^A(3872)) x Br(A(3872) ^ 
TT^ 'K~ J / ip) in the corresponding mode. The number of 
events in the background-subtracted combined distribu- 
tion is 



196.0ti|^. 



For the decay A(3872) Tr+Tr^Tr" J/tA Babar reports Q 

B^^ = [0.6 ± 0.2(stat) ± O.l(syst)] x 10'^, 
= [0.6 ± 0.3(stat) ± O.l(syst)] x 10"^ 

for the charged mode B^ J/ipujK^ and for the neu- 
tral mode B^ — >■ J/tpojK^, respectively. Similarly to 
the two-pion case above, iJa^r stands for the product 
branching fraction Br(B KX{3872)) x Br(A(3872) 
tt+tt^tt'^ J/t/i) in the corresponding mode. The number 
of events in the combined distribution is 



N: 



sig+bg 
3tv 



34.0 ±6.6, N. 



bg 



3tt 



8.9 ± 1.0, 



(2) 



and we assume a flat background. Note, the spectrum re- 
ported in Q and used below appears not to be efficiency 
corrected. However, since only the shape of this spec- 
trum plays a role for the analysis (see number-of-event 
distributions ([TTT) below) and we can reproduce the the- 
oretical spectra of Q, which have the efficiency of the 
detector convoluted in via a Monte Carlo simulation, the 
invariant mass dependence of the efficiency corrections is 
expected to be mild and therefore should not affect our 
analysis significantly. 

Thus the updated ratio of branchings reads Q 



Bs^ Br(A(3872) ^ n+TT-n'Kj/tP) 



B27, Br(A(3872) ^ n+n-J/iP) 
In our analysis we use the ratio ([3]) as well as 
^2^ = 196, A^3^ = 25.1, 



= 0.8 ±0.3. (3) 



(4) 



and he corresponding bin sizes are A£^2ir = 20 MeV and 
AEs^ = 7.4 MeV. 
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III. THEORETICAL tt+tt" 
INVARIANT MASS DISTRIBUTIONS 

As in previous analyses, we assume that the two-pion 
final state is mediated by the p in the intermediate state, 
while the three-pion final state is mediated by the w. It 
was shown in Q that the description of the tt+tt" J/^/; 
spectrum with the 2 ^ assumption is improved dras- 
tically if the isospin-violating p-u mixing is taken into 
account. Theoretical issues of the p-ui mixing are dis- 
cussed, for example, in ,10j . Here we include this effect 
with the help of the prescription used in [ll|, where the 
transition amplitude is described by the real parameter 
A^^p — Aj,^^ — €. Thus, the amplitudes for the decays 
X TT^TT^J/t/j and X tt'^ J / ijj take the form 

+Ax -^j /ipu]Gu)^G pAp-j,2-K , 

+Ax^,j/.,j,pGpeGi^Ai^^^T^, 
where the vector meson propagators are 

Gy^—rni^ — rn^ — irnvrvim)^ V = p,uj, (5) 

with m = TOtttt (m — mTrTnr) being the tt+tt^ (Tr+Tr^vr") 
invariant mass in the tt+tt" J/V' (Tr+Tr'Tr" J/?/;)Jinal state. 
Masses of the vector mesons used below are [12] 



775.5 MeV, 



782.65 MeV. 



The complex mixing amplitude multiplying the uj prop- 
agator used, for example, in [3| to analyse the two-pion 
spectrum, in our notation reads as eG'p(mt^); in partic- 
ular we reproduce naturally the phase of 95° quoted in 
[3]. Note, as we shall only study the invariant mass dis- 
tributions of the two final states, we do not need to keep 
explicitly the vector nature of the intermediate states. 
For the "running" p meson width we use 



p— !-27r 



^ m qinip) 



hp{q{m)) 



hpilimp)) 



where q{m) = ^TO2~~4mf/2, fip{q) = (1 + r^g^) 1/2^ 
with Tp — 1.5 GeV^^ and with the nominal p meson 



width r. 



(0) 



146.2 MeV. 



For the lu meson "running" width (the nominal width 



being H 



(0) 



.49 MeV), the Stt and ir-f decay modes are 



summed, with the branchings 
Br(a; ^ Stt) = 89.1%, Br(cj 
In particular. 



7r7) = 8.28%. (6) 



(7) 



while, for the /l;-i-37r(™)j we resort to the expressions de- 
rived in [isj , with a reduced contact term which provides 
the correct nominal value of the a; — )■ Stt decay width [TJ] ■ 



The transition amplitudes for the decays X — J/ipV 
are parameterised in the standard way, namely. 



Ax~^j/-il>v — 9xvfix{p), 
with the Blatt-Weisskopf "barrier factor" 

foxip) - 1, and hxip) = (1 + r•V)-'/^ 



(8) 
(9) 



for the 1++ and 2 assignment, respectively. Here p 
denotes the J/tjj momentum in the X rest frame. The 
"radius" r is not well understood. If one associates it 
with the size of the X — )■ p{uj)J/'4' vertex, it might be 
related to the range of force. In the quark model this 
radius is r ~ 0.2 fm = 1 GeV~^. This is also in line 
with the inverse mass of the lightest exchange particle 
allowed between J/tjj and p/uj, namely, /o(980). On the 
other hand, a larger value r = 5 GeV~^ is used in the 
experimental analysis of [3|]. Therefore, in the analysis 
presented below we use both values r = 1 GeV~^ as well 
as r = 5 GeV~^, keeping in mind that smaller values of 
r are preferred by phenomenology. 

The theoretical invariant mass distributions for the 
7r"'"7r~ and Tr+Tr^Tr*^ final state take the form: 

dBr27r 



= Bmprp^2rrP^'+' fixip) \RxGp + eGpG. 



dm 



(10) 



^f?xip)\G^ + ^RxG^Gp 



where Rx — gxp/gxu and the parameter B absorbs the 
details of the short-ranged dynamics of the X production. 
The theoretical number-of-event distributions read 



A^27r(m) = X — , 

i>27r dm 

N3n{m) = X 



(11) 



i?37r dm 

The p-uj mixing parameter e is extracted from the 
w — >■ 27r decay width (Br(a; — >• 27r) = 1.53%). The corre- 
sponding amplitude reads 



>2ir — £GpAp^2TT-, 



and we find that 



y^m^mp 



3.4 X 10"^ GeV^ 



(12) 



(13) 



IV. FITTING STRATEGY AND RESULTS 

The number-of-event distributions ([TT]) possess 3 free 
parameters: the "barrier" factor r, the ratio of couplings 
Rx and the overall normalisation parameter B. As out- 
lined above, we perform the analysis for two values of 
r, namely, the preferred value of 1 GeV~^ and a signifi- 
cantly larger value of 5 GeV^^ used in earlier analyses. 
Since the normalisation factor B drops out from the ratio 
of the two branchings, we extract the ratio Rx directly 
from the integrated data, that is from the relation 



c?Br3^ 
dm 



dm 



d&V2T: 

dm 



dm] =S3^/S2^, (14) 
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TABLE I. Sets of parameters and the quality of the combined fits to the data. The last column refers to the presentation of 
the different fits in Fig. \T\ 



Fit 


Wave 




r, GeV-^ 


Rx = gxp/gxu 




CL 


Curve/Markers in Fig. [T] 


S 
Ps 
Pi 


S 
P 
P 


1++ 
2-+ 
2-+ 


5 
1 


'^•^"-'-0.03 

u.uy_o.o2 


1.07 
1.33 
2.77 


37% 
14% 

10-^ 


Solid (red)/Triangles 
Dash-dotted (blue) /Squares 
Dashed (green) /Diamonds 



where the value of the ratio on the right-hand side is 
fixed by Eq. ([3]), and in the integration above we have 
cut off the tt+tt" invariant mass at 400 MeV, as in 3, 
and the tt+tt^tt^ invariant mass at 740 MeV, as in 3- 
Therefore the norm B is our only fitting parameter which 
governs the overall strength of the signal in both chan- 
nels simultaneously, while the shape of the curves is fully 
determined from other sources. 

In TablelH we list the parameters of the 3 combined fits 
to the data, found for the 2 values of the Blatt-Weisskopf 
parameter r. The corresponding line shapes and the re- 
sult of the integration in bins are shown in Fig. [T] 

One can see from Table |T] and Fig. [T] that the best 
overall description of the data for the two channels un- 
der consideration is provided by the 5- wave fit. The 
P-wave fit is capable to provide the description of the 
data of a comparable (however somewhat lower) quality, 
only for large values of the Blatt-Weisskopf parameter r, 
r = 5 GeV~^. The P-wave fit becomes poorer when the 
Blatt-Weisskopf parameter is decreased, and for values 
of r of order 1 GeV~^, the quality of the P-wave fit is 
unsatisfactory, which is the result of a very poor descrip- 
tion of the two-pion spectrum — see the dashed (green) 
curve in Fig.[TJ Varying the ratio of branchings Bj,Tr/B2TT 
around its central value within the experimental uncer- 
tainty interval [see Eq. leads only to minor changes 
in the fits and does not affect the conclusions. 



V. DISCUSSION 

Since no charged partners of the X(3872) are observed 
experimentally, it is supposed to be (predominantly) an 
isoscalar. Then the ratio Rx = gxp/gxui measures the 
strength of the isospin violation in the X — )■ VJ/ip de- 
cay vertex. As discussed above, this ratio is extracted 
directly from the data on the ratio of the branchings ([3]) . 

An isospin- violating observable for a compact charmo- 
nium is the ratio of the branching fractions for the ipi'iS) 
decays into iiJ/t/j and n^J/ijj final states as 



be significant meson- loop contributions [T^ , the estimate 
(|15p is to be regarded as a conservative upper bound for 
the isospin violation strength for compact charmonia. 

In contrast to this, in the S'-wave molecular picture for 
the A", isospin violation is enhanced significantly com- 
pared to the strength (jlSp for it proceeds via intermediate 
DD* states and is therefore driven by the mass difference 
A « 8 MeV of the neutral and charged DD* threshold — 
see, for example, [H, [13] ■ An order-of-magnitude esti- 
mate is provided by the expression 



R 



X 



h{Mx)~Ic{Mx) 



h{Mx) + Ic{Mx) 



0.13, (16) 



where mo is the D meson mass, while Io{Mx) and 
Ic{Mx) denote the amplitudes corresponding to loop 
diagrams with neutral and charged DD* intermediate 
states, respectively, evaluated at the X mass. They are 
composed of two terms, the strongly channel-dependent 
analytic continuation of the unitarity cut, proportional 
to the typical momentum of the meson pair, and the 
weakly channel-dependent principle value term, whose 
size is identified with the inverse range of forces of order 
of 1 GeV (see above). This estimate is within a factor 
of 2 consistent with the value Rx ~ 0.26 found from our 
5- wave fit — see Table HI 

On the other hand, if the X has the quantum numbers 
2 ^, one should expect the isospin violation in the X 
wave function to be of the natural charmonium size, and 
thus of the order of 10^^ — see discussion below Eq. (|15l) . 
since the DD* loop effects are suppressed by the ad- 
ditional centrifugal barrier: the estimate analogous to 
Eq. (dH) now reads Pf"' - (VMdA//3)3 - 2 x lO'^. 
Thus, for the state with the quantum numbers 2 '", one 
expects values of at most Rx ^ 10~^, that are signifi- 
cantly smaller than those following from the data (see Ta- 
ble H]). One is led to conclude therefore that for Rx 0.1, 
needed for the quantum numbers 2 ^ to be consistent 
with the data on the X decays, a new, yet unknown, 
isospin violation mechanism would have to be invoked. 



R 



g^oj,^ _ /Br(V'(25) ^oJ/^) fc3 



g„j/i, V Br(i/^(25) ^ T^J/iP) kl. 



0.03, 
(15) 

where k^^o and k^i are the center-of-mass momenta of tt" 
and 77, respectively, and the tp{2S) branching fractions are 
taken from [T2j . However, since here also the denomina- 
tor violates a symmetry, namely, SU(3), and there might 



VI. CONCLUSIONS 

We conclude therefore that, although the present qual- 
ity of the data in the X — >■ Tr~^ Tr~ tt'^ J / tp channel is not 
sufficient to draw a definite conclusion concerning the 
quantum numbers of the A (3872), the combined anal- 
ysis of the existing two- and three-pion spectra favours 
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FIG. 1. Comparison of the theoretical predictions with the data. Experimental data, taken from 0,01 for the two-pion (left 
plot) and three-pion (right plot) spectrum, respectively, are given by black circles with error bars, while the theoretical results 
are shown both as the continuous number-of-event distributions as well as by markers for the corresponding integrals over 
the bins. Fit 5* is shown as the (red) solid line and triangles, fit Pi is shown as the (green) dashed line and diamonds, while fit 
P5 is shown as the (blue) dash-dotted line and squares. See Table [l] for the details. 



the S-wa.ve fit, related to the 1++ assignment for the 
X(3872), over the P-wave fit, related to the 2 assign- 
ment. We notice that an acceptable P-wave fit calls for a 
large range parameter in the Blatt-Weisskopf form factor 
which meets certain difliculties with its phenomenolog- 
ical interpretation. In addition, while the value Rx — 
gxp/gxu) can be understood theoretically for the 1++ as- 
signment, the value extracted for the 2 ^ assignment is 
too large to be explained from known mechanisms of the 
isospin violation. 
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